mod=10**9+7
t=int(input())
for i in range(t):
  a,b,c,n=map(int,input().split())
  n=pow(2,n,mod-1)
  a=pow(a,n,mod)
  b=pow(b,n,mod)
  c=pow(c,n,mod)
  print((a*b*c)%mod)